Bosonic spectral density of epitaxial thin- film Lai 83Sro.i7Cu04 superconductors from 

infrared conductivity measurements 
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We use optical spectroscopy to investigate the excitations responsible for the structure in the 
optical self-energy of thin epitaxial films of Lai,g3Sro.i7Cu04. Using Eliashberg formalism to invert 
the optical spectra we extract the electron-boson spectral function and find that at low temper- 
ature it has a two component structure closely matching the spin excitation spectrum recently 
measured by magnetic neutron scattering. We contrast the temperature evolution of the spectral 
density and the two-peak behavior in La2-a;Sr2:Cu04 with another high temperature superconduc- 
tor Bi2Sr2CaCu2 0g+5. The bosonic spectral functions of the two materials account for the low Tc 
of LSCO as compared to Bi-2212. 

PACS numbers: 74.25.Gz, 74.62. Dh, 74.72. Hs 



Despite twenty years of extensive research on the high 
critical temperature superconducting cuprates, there is 
as yet no consensus on the pairing mechanism or why 
Tc is so high. It is not known, for example, why 
most of the cuprate materials have a high of over 
90 K while another group, for example La2-xSr2:Cu04 
(LSCO) have a much lower Tc of about 30 K. Many the- 
oretical ideas have been put forward such as the spin- 
charge separation^, preformed pairs 2], exchange of s pin 



fiuctuations[3, 4, 5, 6, 7, 8, .9], and phonons[10. 1 111. 1 12 1 
to list a few of the more prominent proposals. It took 
more than 45 years after the discovery of superconduc- 
tivity in Hg before BCS theory was formulated. Yet less 
than ten years later, accurate and detailed measurements 
of the electron-phonon spectral density a'^F{ijj) responsi- 
ble for the pairing were available from tunneling experi- 
ments for many of the conventional materials. This criti- 
cal function followed from a numerical inversion, centered 
on the Eliashberg equations, of the current voltage char- 
acteristics of tunnel junctions[l3|- Subsequently optics 
was also used successfully to get comparable data[l3, lisj l 
and a'^F{uj) was calculated from first principles and used 
to calculate material specific superconducting properties 
which are often quite distinct from the universal laws of 
the BCS theory 

Many attempts have been made to obtain equiva- 
lent information in the cuprates. An important issue 
in such work is whether or not the idea of boson ex- 
change mechanism and Eliashberg theory can in fact 
be used in highly correlated swtems. Tunneling^, [T3l . 
angle-resolved photoemissionjlll ITtI . [isj . and optics[l9l 
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|22| . |23| . |24| . |25| | have all been used to determine 
an electron-boson spectral density in the oxide supercon- 
ductors but with mixed results. Some argue that the 
resulting spectrum is characteristic of spin fluctuations 



while others attribute it to phonons. Accurate inelas- 
tic neutron scattering measurements of the bosonic spec- 
tra of both phonons and spin fluctuations followed by 
high resolution optical or tunneling spectroscopy on the 
same material would go a long way towards settling this 
issue[2^. A related technique, angle resolved photoe- 
mission (ARPES), is momentum resolved and has shown 
that the boson structure can depend on momentum [27| . 
Thus any definitive comparison with optics would require 
an average of the ARPES data taken at all momenta. 

Very recently inelastic neutron scattering data on the 
local i.e. momentum averaged spin susceptibility in 
LSCO have become available[28] and show two distinct 
energy scales. One peak is centered around 18 meV and 
the another near 40-70 meV with small features extend- 
ing to 150 meV. In this letter we invert datafSi*! on the op- 
tical conductivity of Lai.83Sro.i7Cu04 within an Eliash- 
berg formalism with a Kubo formula for the conductiv- 
ity to obtain the electron-boson spectral density [23|. We 
find a remarkable similarity with the neutron data which 
constitutes strong evidence for the spin fluctuation mech- 
anism. We also compare our new results with previous 
results for BiaSraCaCuaOg+a (Bi-2212). 

For optical spectroscopy La2-a;Sra;Cu04 material 
presents special problems because of the presence of c- 
axis longitudinal phonons in ab-plane optical spectra pos- 
sibly mixed in by an unknown lattice defect. These effects 
are seen in most single crystal spectrajs^l- Fortunately, 
in epitaxial films, these defects seem to be completely 
absent. We have therefore used the low noise data from 
Gao et al. for our refiectance data[2§|. The films are 
820 nm thick, grown on SrTi03 substrate by magnetron 
sputtering. The films show a superconducting transition 
at 31 K with a transition width of 1.5 K. The details of 
the measurements are given in Ref. The raw a {to) 
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FIG. 1: (color online) (Top panel) The optical scattering rate 
l/r°^(tj) in meV as a function of uj (in meV) for five tem- 
peratures. The heavy lines are the experimental data (TO 
phonons have been subtracted to isolate the electronic contri- 
bution). The light solid lines are the fits to the data. (Middle 
panel) The electron-boson spectral density /^x('^) obtained 
from our Eliashberg inversion from the data in the top panel. 
The inset shows the data of VignoUe et al. for a closely re- 
lated sample of La2-i,Sra;Cu04 with x = 0.16 and Tc = 38.5 
K. The data with solid points is at 12 K and a comparison with 
300 K data below 40 meV is also given (open points). (Bot- 
tom panel) The real part of the optical self-energy EJ''(T, tj) 
for LSCO. The arrows show the positions of the sharp peaks 
found in the spectral density at low temperature, w = 15 and 
Lj — 44 meV. Note the sharp rise in the self-energy at these 
frequencies. In the inset we show a simulated self-energy using 
a mode with two Einstein modes. 



data for Lai.83Sro.i7Cu04 [29| shows direct absorption 
by transverse optical phonons as well as the electronic 
background of interest here. After an appropriate sub- 
traction of the phonons, the real and imaginary parts of 
the temperature (T) and frequency (w) dependent elec- 
tronic conductivity a{T,uj) are obtained. For correlated 
electrons, a generalized Drude form applies which defines 
an optical self-energy E°p(T, through the equation 



a{T,Lu) 
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(1) 



tering rate 1/t°p{T,uj) and the optical mass are re- 
lated to E°P(r,tj) by 1/T°f(r,w) = -2E°P{T,uj) and 
uj[m*°P{T,uj)/m-l] = -2I]7(T, w) where 1 and 2 denote 
real and imaginary parts. 

Our method for determining the bosonic spectral den- 
sity begins with the application of the maximum entropy 
inversion 2J| of a simplified convolution integral which 



in UJ- 2E°P(r,w) 
where flp is the plasma frequency. The optical scat 



relates the measured scattering rate 1/t°p{T,uj) to the 
electron-boson spectral density known to be remarkably 
accurate in the normal state. This gives a first numer- 
ical model for /^x(a;) which we then further refine 25| 
through a least squares fit to the optical scattering rates 
using the full solutions of the Eliashberg equations and 
associated Kubo formula. This gives our final model. 
Details about the maximum entropy inversion and the d- 
wave Eliashberg equations can be found in reference [23j| 
and [111. 

In the top frame of Fig. 1 we present results for 
1/t°p{T,lu) vs. u) for five temperatures. All data are 
in the normal state except for the last one (30 K) which 
is just below 31 K, the superconducting Tc of this film. 
The heavy curves correspond to the data which is used 
as the input in the inversion process. Our final results 
for the electron-boson spectral density I'^xi'^) are shown 
in the middle panel. Note the strong temperature de- 
pendence of the low energy structure. When this /^x(cj) 
function is used in the Eliashberg equations, the optical 
scattering rates that result are shown in the top frame 
as the light lines. In all cases the fit to the data is good. 
For the 30 K data in the superconducting state d-wave 
symmetry was assumed for the gap channel. 

The middle frame of Fig. 1 shows the extracted 
bosonic spectral density. Note that at 30 and 50 K both 
spectra show a two-peak structure with peaks at 15 and 
44 meV. As the temperature is increased the spectral 
density evolves towards a single peak which broadens and 
moves to higher energies. We also note tails extending 
beyond 100 meV but these have a low amplitude as com- 
pared to the peaks. All these features are in agreement 
with the data on the local spin susceptibility (shown in 
the inset) measured by polarized neutron scattering in a 
closely related sample by VignoUe et al. [2^ . The data at 
12 K shows the same two peak structure as we have found 
at 30 K. While VignoUe et al. have limited data at 300 K 
they are also in accord with the temperature dependence 
of our spectra: the low energy peak in I^xi'^)^ which ex- 
ists only at low temperature is absent above 100 K as it 
is in the neutron spectra at 300 K. The detailed agree- 
ment between I'^xi^) a-^d the neutron data on the local 
magnetic susceptibility x"{^) is strong evidence that the 
charge carriers in Lai.83Sro.i7Cu04 are coupled through 
a spin Uuctuation mechanism. In related ARPES work on 
highly underdoped LSCO, Zhou et aL[lli] interpret their 
results in terms of phonons. We note that the magnetic 
resonance peak is known to be very weak in highly under- 
doped materials[31] and that (as discussed above) nodal 
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direction ARPES cannot be directly compared with mo- 
mentum averaged optics. 

A phenomenological model for the oxides where spin 
fluctuations replace the usual phonon exchange is the 
nearly antiferromagnetic fermi liquid model (NAFFL) . 
This model is well developed and is anchored in the gen- 
eralized Eliashberg equations. While our approach to 
the inversion of the optical data might not apply in some 
highly correlated metals, it is fully justified in the NAFFL 
model. In particular for a discussion of the applicability 
of Migdal's theorem (a precondition in our approach) we 
refer the reader to the review by Chubukov et al. [3| • 

Two peak structure in the fluctuation spectrum of 
Lai.83Sro.i7Cu04 at low temperature and the disappear- 
ance of the lower peak at high temperatures can be seen 
directly in the optical self-energy itself which follows from 
the raw data without any appeal to microscopic mod- 
els. This can be seen clearly in the bottom panel of 
Fig. 1 where we show our data for — 2I]°^(T, w) for three 
temperatures, T = 30 K in the superconducting state 
(just below Tc = 31 K), at 50 K, above T^, and 250 
K. The heavy black arrows indicate the position of the 
two peaks in I'^xi^)- They coincide with sharp rises in 
—2T,°^{T,oj). Such sharp rises are expected in the nor- 
mal state for coupling of electrons to two Einstein os- 
cillators as shown for a model calculation in the inset 
where a model spectrum J^x(w) has two peaks at 31 and 
62 meV of width 1.2 meV with the second peak having 
twice the spectral weight of the first. The two peaks in 
/^x(ti;) can also be seen directly in the second derivative 
2'!tW{uj) = d'^[(jj/T''P{uj)]/d(jj'^ (not shown here) which is 
a model independent method and is known to be closely 
related to the spectral density in the peak region [l^ . The 
red dashed-dot curve at T = 250 K shows a single rise 
consistent with a single peak in /^x(w) at this tempera- 
ture. 

It is interesting to compare LSCO with similar results 
for optimally doped Bi-2212[25|| (T^ = 96 K). These Bi- 
2212 data are shown in Fig. 2. The optical scattering 
rates on which our inversions are based are shown in the 
top frame while the lower frame shows our results for the 
spectral densities I'^xi^)- Only two temperatures are 
shown 72 K and 300 K for Bi-2212 (dashed curves) and 
30 and 250 K for LSCO (sohd curves). The scattering 
rates are very different in the two materials. The dashed 
curve for Bi-2212 shows a very steep rise with the mid 
point at 92 meV which corresponds to the energy as the 
peak in the real part of the optical self-energy. By con- 
trast the rise in the thick solid blue curve for LSCO is 
more gradual and proceeds on a broader energy scale. 
The rise starts close to zero energy because of the low 
energy peak in J^x(ti>) (thick solid blue curve in the bot- 
tom panel) and is much broader in energy because of 
the second peak around 50 meV. At higher temperatures 
the thin solid and the thin dashed red curves still show 
distinct variations with lu. For Bi-2122 the underlying 
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FIG. 2: (color online) (Top panel) Scattering rates for Bi- 
2212 (dashed lines) and LSCO (solid lines). The two materi- 
als show striking differences at low energies. Scattering sets 
in at very low energies in the LSCO material whereas Bi-2212 
has a negligible scattering below 50 meV. (Lower panel) The 
bosonic spectral function I^xi'^) of the same materials shown 
in the top panel. At low temperature LSCO has two peaks 
while Bi-2212 shows a single peak at much higher energy. As 
the temperature is raised the Bi-2212 peak is weakened dra- 
matically and moves to high frequency ^2^] whereas in LSCO 
the lower peak vanishes while the upper peak remains up to 
250 K. This difference at high temperature can already be 
seen in the top panel: the LSCO sample shows an onset of 
scattering in the 50 meV region while the Bi-2212 curve is 
quite featureless at 300 K. 

spectrum /^x(^) evolved to a single broad peak at 
about 100 meV plus a background. For LSCO the spec- 
trum shows a sharper peak at about 55 meV and the 
background in comparison is small. This translates into 
a room temperature scattering rate in Bi-2212 which is 
featureless and flat while in LSCO the bosonic peak is 
very much present giving rise to a marked shoulder in 
the scattering rate in the 50 meV region. 

What are the implications of these spectra for super- 
conductivity? Tc in Bi-2212 is three times larger than 
in LSCO. For an s-wave electron-phonon superconductor 
the same spectral density /^x(^) enters both renormal- 
ization and gap channels of the Eliashberg equations. In 
a d-wave superconductor, however, the spectral density 
that enters the gap channel is a d-wave projection of the 
electron-boson exchange process rather than its s-wave 
projection which we have determined from optics and 
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shown in the middle frame of Fig. 1 and the lower frame 
of Fig. 2. For simplicity we can assume that these two 
quantities differ mainly by a numerical factor 33| • This 
factor can be fit to the known value of the critical temper- 
ature Tc which we determine from the complete numeri- 
cal solutions of the d-wave Eliashberg equations. We find 
that Tc is well represented by the simple modification of 
the McMillan equation [H Hi] 



for preparing the film. 



ksTc'^h LOin exp 



(2) 



where 1 -f A* comes from the normal state renormaliza- 
tion of the dispersion curves and A'' is the interaction in 
the gap channel. Here u;;„ is the average boson energy 
defined by Allen and Dynes [35] and is the same for s- or 
d-channel cases. 

Using our functions we find that the A"* values 

are nearly the same for Bi-2212 (1.85) and LSCO (1.90). 
However the value of a;/„ differ by a factor of two, '^50 
meV for Bi-2212 and -25 meV for LSCO. Therefore the 
softening of the spin fluctuation spectrum in LSCO as 
compared to Bi-2212 accounts for a factor of 2 difference 
in Tc- The remaining difference is traced to the value of 
A* which is larger in LSCO (3.40) as compared to Bi-2212 
(2.50). The renormalization factor 1-1- A* in the modified 
McMillan equation (Eq. (2)) is pair breaking and this 
accounts for the rest of the difference in Tc values. While 
we do not predict Tc we can explain in a robust way 
from our analysis the factor of 3 difference in Tc between 
LSCO and Bi-2212 samples. 

In summary we have measured the fluctuation spec- 
trum in LSCO, I'^xi'^) various temperatures. At low 
temperature it shows two characteristic energy scales in 
remarkable agreement with the local (q averaged) spin 
susceptibility recently found in polarized inelastic neu- 
tron scattering experiments. As the temperature is in- 
creased, the low energy peak disappears in accord with 
the neutron results. While the maximum entropy tech- 
nique and least squares fit to the measured optical scat- 
tering rates within an Eliashberg framework is employed 
here, the two-peak structure in /^x(a;) can be seen di- 
rectly in the raw data for the real part of the optical 
self-energy. In contrast optimally doped Bi-2212 reveals 
a very different behavior showing a single sharp peak 
centered at 60 meV with a valley above it and a broad 
low intensity background extending to energies up to 400 
meV. Finally, the bosonic spectra derived from our anal- 
ysis, fully account for the low superconducting transition 
temperature of LSCO as compared to Bi-2212. 

Acknowledgments 

This work has been supported by the Natural Science 
and Engineering Research Council of Canada and the 
Canadian Institute for Advanced Research. FG and DBT 
acknowledge support from the NSF and DOE through 
grants DMR-0305043 and DE-AI02-03ER46070." Also 
we want to acknowldege J. Talvacchio and M.G. Forrester 



[2] 
[3] 

[4] 

[5] 

[6] 

[7] 
[8] 

[9 

[lo; 
[11 

[12 
[13 

[14 

[15 

[16 
[17 
[18 
[19 
[20 

[21 
[22 
[2Z 

[24; 

[25; 

[26' 
[27 
[28' 
[29 
[30 

[31 

[32; 



[33 
[34; 

[35] 



Electronic address: hwangjs@mcma ster.ca| 
P. W. Anderson, The theory of superconductivity in the 
high-Tc cuprates. Princeton University Press, Princeton, 
New Jersey (1997). 

V. J. Emery and S. A. Kivelson, Nature (London) 374, 
434 (1995). 

A. V. Chubukov et at, The Physics of Conventional and 
Unconventional Superconductors edited by K.H. Benne- 
mann and J.B. Ketterson, Springer- Verlag (2003). 
J. P. Carbotte, E. Schachinger and D. N. Basov, Nature 
(London) 401, 354 (1999). 

Ar. Abanov and A. V. Chubukov, Phys. Rev. Lett. 83, 
1652 (1999). 

J. F. Zasadzinski et ai, Phys. Rev. Lett. 96, 017004 
(2006). 

J. C. Campuzano et al, Phys. Rev. Lett. 83, 3709 (1999). 
M. R. Norman and H. Ding, Phys. Rev. B 57, R11089 
(1998). 

P. D. Johnson et al, Phys. Rev. Lett. 87, 177007 (2001). 

A. Lanzara et al, Nature (London) 412, 510 (2001). 
X. J. Zhou et al, Phys. Rev. Lett. 95, 117001 (2005). 
J. Lee et al. Nature (London) 442, 546 (2006). 

W. L. McMillan and J. M. Rowell, Phys. Rev. Lett. 14, 
108 (1966). 

B. Farnworth and T. Timusk, Phys. Rev. B 14, 2799 
(1974). 

F. Marsiglio, T. Startseva and J. P. Carbotte, Phys. Lett. 
A 245, 172 (1998). 

J. P. Carbotte, Rev. Mod. Phys. 62, 1027 (1990). 

A. A. Kordyuk et al, Phys. Rev. Lett. 97, 017002 (2006). 

T. Valla et al, cond- mat/0610249 (2006). 

R. T. Collins et al, Phys. Rev B 39, 6571 (1989). 

A. V. Puchkov, D.N. Basov and T. Timusk, J. Physics, 
Condensed Matter 8, 10049 (1996). 

J. J. Tu et al, Phys. Rev. B 66, 144514 (2002). 

S. V. Dordevic et al, Phys. Rev. B 71, 104529 (2005). 

J. P. Carbotte and E. Schachinger, Ann. Phys. 15, 585 

(2006). 

E. Schachinger, D. Neuber and J. P. Carbotte, Phys. Rev. 
B 73, 184507 (2006). 

J. Hwang, T. Timusk, E. Schachinger and J. P. Carbotte, 

Phys. Rev. B 75, 144508 (2007). 

M. Chiao, Nature Physics 3, 148 (2007). 

T. Sato et al, J. Phys. Soc. Japan, 76, 103707 (2007). 

B. Vignolle et al. Nature Physics 3, 163 (2007). 

F. Gao et al, Phys. Rev. B 47, 1036 (1993). 

M. Reedyk and T. Timusk, Phys. Rev. Lett. 69, 2705 
(1992). 

Z. Yamani et al cond-mat/0706.0327. 
J. P. Carbotte and E. Schachinger, Model and Methods 
of high-Tc Superconductivity Edited by J. K. Srivastava 
and S. M. Rao, Nova Science Publishers Inc., Hauppauge 
N.Y. p 73 (2003). 

W. L. McMillan, Phys. Rev. 167, 331 (1968). 

P. J. Williams and J. P. Carbotte, Phys. Rev. B 39, 2180 

(1989). 

P. B. Allen and R. C. Dynes, Phys. Rev. B 12, 905 (1975). 



